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Two utmost cases of super-extreme event's influence on the velocity autocorrelation function 
(VAF) were considered. The VAF itself was derived within the hierarchical Weierstrass-Mandelbrot 
Continuous- Time Random Walk (WM-CTRW) formalism, which is able to cover a broad spectrum 
of continuous-time random walks. Firstly, we studied a super-extreme event in a form of a sustained 
drift, whose duration time is much longer than that of any other event. Secondly, we considered 
a super-extreme event in the form of a shock with the size and velocity much larger than those 
corresponding to any other event. We found that the appearance of these super-extreme events 
substantially changes the results determined by extreme events (the so called "black swans") that 
are endogenous to the WM-CTRW process. For example, changes of the VAF in the latter case 
are in the form of some instability and distinctly differ from those caused in the former case. In 
each case these changes are quite different compared to the situation without super-extreme events 
suggesting the possibility to detect them in natural system if they occur. 

PACS numbers: 89.65.Gh, 02.50.Ey, 02.50.Ga, 05.40.Fb, 02.30.Mv 



I. INTRODUCTION 

One of the most remarkably emerging observation within the natural and socio-economical sciences is that the 
empirical data which they supply are frequently punctuated by rare extreme events or black swans which can play a 
dominant role. This observation is usually quantified by power-laws or heavy-tailed probability distributions of event 
sizes (cf. pHIII and references therein). However, as it was pointed out in [TjJ, there is an empirical evidence that 
something important beyond power-laws does exist. In this context, the concept of super- extreme events, outliers or 
dragon-kings was introduced. 

By the super-extreme event, outlier or dragon-king 1 we mean an event of the values that are abnormally different 
than values of other events in a random sample taken from a given population [T^J. It is therefore an anomaly, an 
event usually removed in order to obtain reliable statistical estimations. The term "outlier" emphasizes the spurious 
nature of these anomalous events, suggesting discarding them as errors or misleading monsters. In contrast, the 
term "dragon- king" emphasizes their relevance in the dynamics and the importance of keeping them to be able to 
understand the generating process. The super-extreme events are statistically complementary to extreme events, as 
documented in [13|. Moreover, the idea that dragon- kings are often associated with the occurence of a catastrophe, 
a phase transition, and bifurcation as well as with a tipping point, whose emergent organization produces visible 
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1 The poetic term "dragon-king" stresses that: (i) we deal with exceptional event which is a completely different kind of "animal" in 
comparison with the usual ones concerning the rest of events, (ii) it emphasizes the importance of this event being outside the power-law 
like an absolute monarchy standing above the law or as the wealth of a monarch owning a finite part of the whole country, beyond the 
Pareto distribution of its citizen's wealths 1 14|| . 
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precursors was also developed there. 

The main goal of the present work is to analytically demonstrate and numerically simulate two utmost cases of 
the influence of dragon-kings on the velocity autocorrelation function (VAF) of a random walker. Herein, we studied 
VAF in the frame of the Weierstrass-Mandelbrot Continuous-Time Random Walk (WM-CTRW) formalism developed 
in [R||la | . This model is a hierarchical version of the canonical Continuous-Time Random Walk (CTRW) formalism 
[19l - l2l| in which the hierarchical spatio-temporal waiting-time distribution (WTD) was assumed as its basic quantity 
(cf. Section fiTI Al in this work as well as Equation (20) in (l7l|). 

Note that the WM-CTRW formalism is sufficiently generic and flexible. It is able to cover various types of diffusion, 
i.e. from normal diffusion, through the superdiffusion (e.g. the persistent fractional Brownian motion, fBm, (22.]), and 
the ballistic one, to the Levy walk. Remarkably, the extreme events are contained in the spatio-temporal structure 
of time series obtained within the WM-CTRW formalism by stochastic simulation. This is explained in details in 
Section |lTJ The extreme events moderate the relaxation of the system to equilibrium (or partial equilibrium), e.g. by 
changing the relaxation from exponential to power-law. Furthermore, we used the WM-CTRW formalism because 
we (superficially) verified that the hierarchical WTD described quite well its empirical counterparts obtained for 
continuous quotation of the exchange rates on a currency exchange market as well as of the share price trading on a 
stock exchange. 

The present paper is organized in the following manner. First (in Section [TTJ) , our problem is defined along with 
definitions of the most relevant quantities used in our study. Then (Section the WM-CTRW formalism is defined 
and, next, a superdiffusion phase is discussed. Subsequently, Section IIVI presents derivation of VAF including a 
dragon-king event. The comparison of predictions of our theoretical formulae for the VAF with the results of the 
simulations is contained in Section |VJ Finally, summary and concluding remarks are presented in Section IVII 



II. DEFINITION OF THE PROBLEM 



In the present work we consider, as an example, a superdiffusion case within the WM-CTRW formalism, where 
so-called weak ergodicity is obeyed [23|-[25[; that is, the mean value of the waiting-times between turning points of a 
random walk trajectory is finite. As the probability of appearance of a dragon-king is extremely small (that is, the 
waiting time for its appearance is too long for all practical purposes), it was produced after the simulation and next 
put "manually" somewhere inside the series. In that sense, the dragon-king's appearance can be considered as an 
exogenous event. 

We answer the question of how much the stationary velocity autocorrelation function, C(At), derived from a given 
time series is changed when this time series is suddenly punctuated by a single-step super-extreme event. We consider 
two types of super-extreme events: 

(i) the long-drawn event which has super-extremely long duration time td (cf. Figure [T]) and 

(ii) the shock, or sudden jump, of a random variable X, which has super-extreme velocity Vd (cf. Figure [2]). 

The background of the corresponding dragon-kings' definition is again presented in Figures [T] and [H The resulting 
VAF, involving a dragon-king, is denoted below by Cd(At). 
The VAF is defined as 

VAF(Ai) = (v(t r ) v(t' + At)) - (v(t')) (v(t' + At)} 

(v v ) (At) (v ) (v ) — < * n ^ e a b sence °f a dragon- king, 

— [viv 2 ){ ) \vi) [v 2 ) — I Cd(At), in the presence of a dragon-king, ^ ' 

where (. . .) means the moving average (or averaging over current time t' < t to t — At) in a given time-window At (cf. 

Figure [T] and [2] and considerations in Section ITV|) and v(t') d =' [X(t') — X(t' — dt)]/dt, where the time-discretization 
step is dt < At, t' . The velocities v\ and v 2 are defined here at the beginning and at the end of the time- window At, 
respectively. 

Both quantities, C(At) and Cd(At), are studied analytically (in Sections Mil and HV|) and by numerical simulations 
(in Section |V]) because our task is to find relations between C(At) and Cd(At) and verify them by simulations for 
cases (i) and (ii) mentioned above. 



III. WEIERSTRASS-MANDELBROT CONTINUOUS-TIME RANDOM WALK 



The task of this section is to briefly sketch elements of the WM-CTRW formalism that are useful for our analysis. 
An explicit form of the VAF, C(At), in the absence of a dragon-king was described in details in our earlier works 
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FIG. 1: Schematic time series for the total time ttot containing the dragon-king event (represented by the longest sloping straight 
line) having constant velocity Vd and duration time td- The time series before and after the occurrence of the dragon-king have 
duration times II and in, respectively. Obviously, t to t = £l + t d + tR. 



Definition of the WM-CTRW formalism 



The Weierstrass-Mandelbrot Continuous-Time Random Walk is defined by the hierarchical spatio-temporal waiting- 
time distribution (WTD). This WTD 2 is given by the following weighted hierarchical series 

oo 

^(x,t)=Y / w(j)^(x,t), (2) 

3=0 

where x is the walker single-step spatial displacement passed (with constant velocity) within the time interval t and 
the weight w(j) is given by the probability distribution 



W (j) = ^j(l-^),N>1, ,=0,1,2,.... (3) 

This weight can be interpreted as the probability of exactly j consecutive successes in some Bernoulli series, where 
l/N is the probability of a single success and the conditional single- level WTD is assumed in the factorized form of 
two different single- variable distributions / and h 

M*, *) = -^f( — ) —h(-). (4) 



2 The complete definition of the WM-CTRW model additionally requires a special treatment of the first step of the random walker |17| . 
However, it is irrelevant when moving-average is performed. Therefore, we do not consider this special treatment in this work. 




FIG. 2: Schematic time series containing a dragon-king in the form of a shock (represented by the longest vertical straight 
line). The time series before and after the dragon- king have duration times tL and tit, respectively while the duration time of 
the shock is t# = dt being so short that it cannot be visualized in the plot (t to t = tt + dt + tit). 



In Expression we used a simple representation of a random walk (but not random jumps or flights). That is, 
conditional temporal and spatial probability distributions we assumed are of the forms 

,(111)^(111^) (5) 

and 



vqvH ) 2 \vqvH 



cxp , (6) 



T T3 ) \ T T3 

respectively 3 . The mean duration time, tqt^ , of the random walker single step and its velocity, vqv , are both associated 
with the j th level of the spatio-temporal hierarchy. This level is the same for temporal and spatial partial probability 
distributions, introducing the spatio-temporal coupling in ([2]). Herein, we marked the calibration parameters by to 
and vq. 

As it is seen from @, our WM-CTRW formalism belongs to the non-separable category of CTRW, which makes 
the diffusion phase diagram quite rich flU Il7|. 

From ([2]), it is easy to derive useful temporal and spatial single-step moments 

/■OO rCG ( X — l/N r -I 

(t)= dtt dx^(x, t)=T \ ^7N' 7 ° > L > (7) 
Jo J-oo oo, for a < 1 



where a = In N / In r and 



<* 2 ) - r dxx* Halt^x, t) = (6 ) 2 { 2 ^W' f L [> 2 ' (8) 
J-oo Jo oo, for 13 < 2 



Obviously, detailed forms of the scaling functions / and h are less important for macroscopic displacement and asymptotic long time, 
respectively. Here, we used their simple explicit forms to make our calculations easier. 
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where b$ = vqTq, b = vt and /3 = In N / In b. 

Notably, a and (3 are parameters which control the different phases of diffusion. It was proved in [TtJ that the (total) 
fractional diffusion exponent 77, which governs the asymptotic time-dependence of the walker multi-step mean-square 
displacement (cf. Equation (11) in jl7j). 

<* 2 (A*)> « (At)", (9) 

itself depends on a and /3 (cf. Table 1 and diffusion phase diagram shown in Figure 1 in [l7|). We term these parameters 
temporal and spatial diffusion exponents 4 , respectively. The quantity D st is the so-called fractional diffusion coefficient 
emphasizing its association with the stationarized random walk [17]]. It depends on partial (temporal and spatial) 
diffusion exponents (see also Table 1 in [ItJ )• Notably, the total time is ttot = an d displacement of the walker is 
X(t tot ) = J2 X (where ^ means the sum over successive steps of the walker). 

Expressions © - @ enable simulation of a random walk trajectory (schematically shown in Figures [1] and ^ in 
continuous time. This is because they define the corresponding stochastic dynamics (considered in Section IV Ap . 

B. Superdiffusion phase 

The WM-CTRW defines a stationary stochastic process valid only for the case where the mean waiting-time, (t), 
is finite, i.e. for the case where the temporal exponent a > 1 (cf. Equation (JT]) as well as Equations (4) and (5) in 



17 1). Under such conditions, the diffusion exponent r\ = 2H, where < H < 1, is the well known Hurst exponent 
22j . All our analytical calculations are confined to the case where the (multi-step) mean-square displacement is finite 
for a finite time and superlinearly increases with time for asymptotically long time. That is, we are confined to the 
superdiffusion phase where 77 > 1. This regime is of interest because dragon-kings could be argued to be least relevant 
in such phase. Our goal is to demonstrate that dragon-kings make a significant impact even in such a superlinear 
phase. 

The superdiffusion phase is restricted by the following inequalities 

i<i<i + £- CM) 

2 p 2 la 

For the superdiffusion phase, we can easily derive the VAF, in the absence of a dragon-king, in the form (cf. 
Equation (13) in 0) 



where the fractional diffusion coefficient is given by 

1 - 
InJV 



1—4? na 



sin 

and the fractional diffusion exponent is 

"= i+2 «(H)- (13) 

As the fractional diffusion exponent 77 < 2, the velocity autocorrelation function, C(At), given by Expression (p~T|) 
vanishes for extremely long At. In Section HVl we prove that the presence of a dragon- king in the time series can lead 
to violation of this property. 

IV. DERIVATION OF FORMULAE FOR C d {At) 

In this section we derive relations of two different types between the estimator of Cd(At) and the estimators of 
some VAFs concerning time series in absence of any dragon- king. That is, in Section lIV Al we consider case (i) while 
in Section [IV Bl case (ii), where both were already defined in Section HT1 



4 It can be also proved that there exists a simple exponent, which is a function of the partial diffusion exponents a and /3, that governs 
the time-dependence of the walker single-step mean-square displacement, ([x(t)] 2 ) = dxx 2 ip(x, t), for asymptotic time interval t. 
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A. C<j(At) for the case of sustained dragon-king 

In this section we assume that 

(1) dt<At< At max < t L , t R and 

(2) At max < td, where At max is a maximal value of At. 

Importantly, assumption (2) is violated for case (ii) defined in Section HTl above and considered below in Section [IV Bl 
We divide the expected value of the product of velocities, (vi V2) (At), present in Formula (fl), into five weighted 
different components. They are estimated by the following terms 

M=5 

( Vl v 2 ) (At) = («i «2> m (At) w m . (14) 

m— 1 

It is straightforward to derive each component (v\ V2) m (At) separately by using the proper moving-average estimator. 
The corresponding weight w rn is easy to obtain. 

(1) The first component (for m = 1) is defined for the case, where both velocities v\ and V2 are placed before the 
dragon-king position. This component, together with the corresponding weight are, as follows 

1 tL ~ At f — At 

<»' «>■ < At > = irhn £ »«' + At >- « * ^fk (15) 

t'—dt 

and relates to a random walk in the absence of any dragon-king. More precisely, this component represents the 
situation where velocity v\ is placed inside the time interval [0, t^ — At] while velocity V2 can be placed both 
inside this time interval or at its right border t^ — At. 

(2) For the second component (m = 2), velocity v% is placed before the position of the dragon-king while velocity V2 
equals the velocity of the dragon-king Vd- The second component and the corresponding weight takes the form 

(viv 2 ) 2 (At) d = i v(t')v d = (v) L v d , W2 = • TTAi ' ( 16 ) 

where Vd is the velocity of the dragon- king (or the slope of the longest straight line shown in Figure [T|) . This 
term describes the first cross situation where velocity V\ is located inside the time interval [t^ — At, t^] or at its 
right border while velocity V2 is placed inside the time-interval [tL,tL + At] or at its right border. Obviously, 
Vd is the constant velocity of the dragon-king. 

(3) This case is particularly simple as both velocities v\ and v 2 are equal to the dragon-king's velocity Vd- The third 
component (m = 3) is the simplest one and together with the corresponding weight assume the forms 

(v lV 2) 3 (At) d ^vl W3 ^-l^- (17) 

Hot — ^ 

and corresponds to the case where both velocities v\ and 1)2 are placed inside the dragon-king's time interval 
[tL,tL + td] or velocity V2 can be also located at its right border. 

(4) Herein, velocity vi equals the dragon-king's velocity Vd while velocity V2 is placed after the position of the 
dragon-king. The fourth component (m = 4) and the weight are 

, , 1 t L +t d . 

( VlV2 ) 4 (At) d ^— v d v(t')=Vd(v) R , W4 d ^-—^-, (18) 

t'=t L +t d -At+dt tot 

corresponding to the second cross case. Precisely, for this case velocity v\ is placed inside the dragon-king 
time interval \pL + td — At, t^ + td] or at its right border while velocity V2 is placed inside the time interval 
[tL + td, tL + td + At] or at its right border. 
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(5) For the fifth component (m = 5), positions of both velocity v\ and v 2 are placed after the position of the 
dragon-king. The fifth component and the weight are, as follows 



(vi v 2 } 5 (At) 



dcf. 



1 



t tot - At 

E 



t R - At 

t'=t to t—tR+dt 



v(t')v(t' + At), w 5 



dcf. t R — At 

~~ At' 



ti 



(19) 



being analogous to those of the first component and weight. More precisely, they are defined for the case where 
both velocities V\ and v 2 are placed inside the time interval [ttot — t R ,ttot] (without of the dragon-king) and 
velocity v 2 can be also counted at time t to t- 

Subsequently, we can calculate the estimator, which already includes the full influence of the sustained dragon-king 
C d (At) = (v 1 v 2 )(At)-(v 1 ) (v 2 ) 

jL-At/ttot, , . 1R- At /ttot i \ . At I ttot ,i \ , / \ \ , 7rf ^ At /t tot 2 

(ui v 2 ) l (At) + — ( Vl v 2 ) 5 (At) + ((%+ (v) R ) vd + xTu — y d 



1 - At /ttot 
1 

" (1 - At/t tot ) 2 
by using definitions 



1 - At/tt 



1 - At/tt 



1 - At/t t 



[lL (v) L + (lR- At /t tot ) (v) R + idVd] [(lL - At /ttot) (v) l + 1h( v )r + 7d«d] 



(20) 



«2 



dcf. 



dcf. 



1 



At 



At 



' t L ttot- At 

J2 v(t')+t dVd + v{ t') 

.t' — dt t' —ttot —tR-\-dt 

tL ttot 

t'=At+dt t'=ttot-t R +dt 



ttot ~ At 
1 



[*l (v) L + (t R - At) (v) R + t d v d ] 



t 



tot 



At 



[(ti - At) (v) L + t R (v) R + t d v d ] 



(21) 



where (v) L and (v) R are partial mean velocities defined by the random walk on the left and right hand side of the 
dragon-king, respectively (cf. Figure[T|). We set velocities (i>i) L and (vi) R equal to zero as no drift is present in the 

system. Besides, we used dimensionless parameters 7^ d =;' ti,/ttati 1r d = t R /t to t and j d =' t d /t to t- 

Because drift is absent in the system, Equation (|2H)) . taking into account Definitions (|2"Tj) . assumes a simpler form 

n ( A+^ lL—At/ttot , v / * .\ . 1R- At /ttot I . 
Cd(At) = irr77 — ( Vl v 2 ) 1 (At) + — — — ( Vl v 2 ) 5 (At) 

1 — LAI I Hot 

Id \ At /ttot 



1 - At /t t ot 
Id 

1- At /hot 



1 



1 - At/tt 



1 - At/t t 



(22) 



Apparently, the quantity C d (At) depends (in general) not only on the relative variable At /ttot but also on the 
parameters jl, j R) and "/ d . Hence, this quantity depends on the position of the dragon-king inside a time series, i.e. 
this quantity is, in general, a non-stationary one. However, it does not depend on the sign of the dragon-king velocity. 
Moreover, the origin of the dragon-king is irrelevant. That is, our derived formula is sufficiently generic in the sense 
that it is valid not only for the WM-CTRW but for any random walk. 

For sufficiently wide time window AtMAX which still obeys AtMAX /ttot "C 1, the quantity C d (At) simplifies into 
the asymptotic formula 



C d (At) 



Id 



1 - At/U 



1 - 



Id 



1 - At/U 



At/U 



1 - At/t t 



(23) 



As long as the dragon-king is present, C d (At) does not vanish, even though estimators (vi v 2 ) 1 (At) and (v\ v 2 ) 5 (At) 
(which can strongly fluctuate) are decaying. 
Obviously, Formula (|23p takes a simpler form 



C d (At)^ ld (l- ld )v 2 d 



if a strong but reasonable inequality 



At 



MAX 



< min(7i, 7_r, j d (1 - jd)) 



(24) 



(25) 



is obeyed. 
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1. Important case of (vi V2) x (At) = (vi 1*2)5 (At) 

By assuming that random walks before and after the dragon-king appearance within the time series are statisti- 
cally identical (although corresponding trajectories could be quite different) we can write, by neglecting unavoidable 
fluctuations, that (vi v 2 ) x (At) = (vi 772)5 (^) = C(Ai). These equalities enable transforming Expression (|2"2"j) into a 
stationary form 



l_^_2AtA^ 
V ' 1-At/ttot 



Id f 1 _ Id \ At/t t 



,2 



I'.i 



(26) 



1 - At/t to t V 1 - At/t tot J 1 - At/t t 
which is our reference formula. This stationary form is our achievement which enables several applications of Formula 



The main difference between C(At) and C^At) is that the former asymptotically vanishes while the latter does 
not. This difference provides a tool which allows for distinguishing a power-law relaxation, controlled by rare extremes 
or black swans, from a decay controlled by the dragon-king. Indeed, we consider Expression (|2"o]) as a reference one, 
also relevant for more complex cases. 

If Assumption (J2HJ) is valid, further simplification of Expression ([2"6")) can be made 

C d (At)^(l- ld )[C(At)+ ld vj}. (27) 

This expression depends on two parameters 7^ and v d fully characterizing the dragon-king. These parameters can be 
easily determined, e.g. from the initial and asymptotic nonvanishing values of C d (At). In fact, predictions of this 
simple formula is compared in Section [V] with corresponding results of simulation. 



B. Cd(At) for the case of a shock (case (ii) of Section [TTJ) 

To consider the case of a shock within the time-series, we replace Assumtion (2), given in Section llV Al by the 
following one 

t d = dt. (28) 

That is, the duration time of the dragon-king is assumed to be as short as possible and equals to the time discretization 
step dt. The Assumption (|2"5|) imposes a modification of components m = 2 and to = 4 in Expression (TT4"]) . Namely, 
components to = 2 and m = 4 together with their corresponding weights are replaced by expressions 



dcf. 



dt 



\viv 2 ) 2 = v L v d , w 1 = - — , 

Hot — L±t 

and 

dcf. dt 



{.ViV 2 ) A = v d v R , w 2 



hot - At 

respectively. Furthermore, instead of component to = 3 we have 

At - dt 

{viv 2 ) 3 , w 3 = - — , (29) 

Hot — L±l 

where the current velocities v\ = vl and v 2 = vr are now located before and after the position of the shock, 
respectively. Hence, a new relation for C d (At) has the form 

C d {At) = [{t L - At) {v x v 2 ) 1 (At) + dtv L v d + (At - dt) ( Vl v 2 ) 3 (At) + dtv d v R + (t R - At) ( Vl v 2 ) 5 (At)] - — l —- 

Hot — 

- [t L (v) L + dtv d + (t R - At) (v) R ] {(t L - At) (v) L + dtv d + t R (v) R \ (30) 

(ttot — Lit) 

By analogy to case (i) of a sustained dragon-king, we set (vi v 2 ) 1 (At) = (vi v 2 ) 3 (At) = (v± v 2 ) 5 (At) = C(At) that, 
in the case of no drift present in the system, simplifies ([3U]) into the form 

Cd(At) = ^- A *- 2d * C(At) + ^ + v»l X * _ (^r) 2 , (31) 
v ' ttot -At K ' Uot -At \tt r ,t-At v ' 



9 



where Xd = dtvd is the value of the shock. Note that vl and vr are velocities of the walker, which are separated 
by the time interval 2At. From definition of the time step dt we have dt <C t to t — At. Hence, Equation (|31[) again 
assumes a simpler form 

a(a() = c(a()+ (^i^_(_ii_) 2 , (32) 

which is our second basic formula. This formula is further modified in Section [V] to make Cd(At) better suited for 
comparison with our results obtained by simulations. 



V. ALGORITHM AND RESULTS 



A. Stochastic dynamics 

The stochastic dynamics or algorithm simulating successive single-step displacements of the walker following a 
WM-CTRW consists of three stages. 

(i) The drawing of the level index j from distribution §3§ in each spatio-temporal step separately. 

(ii) The calculation of the duration time t of the single-step (or its elapsed time) from the stochastic equation 

t = -t t 3 ln(l - R), (33) 

where Re [0, 1[ is a random number drawn from the random number generator (or from the uniform distribution 
confined to a unit interval); note that Equation Q33[) is equivalent to Equation ((SJ) (after an application of well 
known method of the cumulative distribution function inversion). 

(iii) The determination of the single-step displacement by Equation 

x(t) = £v v j t, (34) 

where stochastic variable £ is a dichotomic noise (i.e. £ = +1 or —1 with equal probability 1/2) and xit) d ==' 
X(t') — X(t' — t), where t' is a current time (and not a time interval). 

In this way construction of a single continuous in time trajectory of a random walk is possible. Obviously, the 
parameters tq,t, vq, and v were fixed at the beginning of a whole simulation. Further in the text and in all our 
simulations we set the calibrating parameters tq — 1 and vq = 1. In fact, this trajectory is constructed in the frame 
of a non-stationarized version of a CTRW, where there is no special treatment of the initial step (l(l UtI l20l | . That 
is, the hierarchical spatio-temporal waiting-time distribution {J2j> f° r the initial step is the same as for all other steps. 
Indeed, the stationarized version of the CTRW, which differently treats the waiting-time distribution for the initial 
step, corresponds to the WM-CTRW formalism 5 [l7j. The VAF obtained theoretically within the WM-CTRW enables 
for comparison with the corresponding VAF obtained from simulations. This is possible because in simulations we 
deal with moving-averages, which by definition, average over the initial state, thus supplying the required asymptotic 
stationary VAF. 

The above given algorithm was used in Section IV Bl to simulate the required basic random walk trajectory. The 
trajectory constructed in such a way is punctuated "manually" by the sustained dragon-king or by the shock dragon- 
king. However, preparation of the former dragon-king requires some explanation. 



1. Preparation of sustained dragon-king 



The stochastic dynamics defined by Equations (|33[) and (|34[) is controlled by three random variables: j, R, and £, 
while only the integer level index j, drawn from the distribution ([3]), and random number R are responsible for the 



5 It can be proved that the WM-CTRW formalism, even asymptotically, is more general than the fractional Brownian motion introduced 
by Mandelbrot and van Ness |22|| . It is because the propagator derived within the WM-CTRW formalism can be non-Gaussian (for 
1//9 > 1/a). 
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size of a single step. For the case of the sustained dragon- king this dynamics is simplified by replacing Equation (l33l) . 
which simulates an exponential distribution of interevent times ([5]) , by the corresponding reduced expression 

t = t d = T T 3d . (35) 

That is, the size of the sustained super-extreme event is controlled only by a single random variable j = jd as are the 
velocity of the sustained dragon-king Vd — vq v 3d , its duration time t ( i = tq r Jd , and displacement Xd = £ £>o b 3d , where 
bo = vq tq and b — vt. 

We can say that the stochastic process defined by Equations (|3"S")l and is simplified. This process is a dicrete 
in time and hence in space, where 

x(t) = X d = £ v Q tq v jd T jd = £ b b jd . (36) 

If the index jd of a given dragon-king is fixed, the interevent times for this process cannot fluctuate. Its discrete-time 
step equals to the mean-time defining the exponential distribution (|33[) or ((SJ). 

Apparently, we deal with two stochastic processes: (i) the first one which prepares the WM-CTRW trajectory or 
stochastic spatio-temporal hierarchy of events and (ii) the second process which generates the sustained dragon-king 
from the simplified, discrete in space and time hierarchical random walk. In fact, the latter process is used in Section 
IV A 21 to illustrate the definition of black swan. 

Right now, it is easy to separate in simulation the sustained dragon-king from the rest of events belonging to the 
stochastic spatio-temporal hierarchy of events (cf. Section IV A 2|) . Namely, it is sufficient to choose the level index jd 
much larger than the maximal level index j = juAX of the hierarchy defining the longest single-step displacement 
of any other event. The above described simplicity is the main reason for such a way of selection of a sustained 
dragon-king. 



2. Hierarchical random walk and black swans 



It is decisive for our considerations that the ratio of successive weights 

w(j + 1) = l_ 
w(j) N' 



(37) 



is already j independent. This means that steps defined by the level index j are TV times more likely than those of the 
next step of larger order j + 1. Therefore, one expects (on the average) that the walker will perform N 3 shorter steps 
before performing the next step of larger order. Hence, we can explain how extreme events or black swans control the 
hierarchical spatio-temporal structure of events in the frame of a simplified hierarchical random walk. 
Here we consider, as a typical quantity, the mean-square displacement of the process 




(X 2 )(L) = ( V> )=L(x 2 ), (38) 



where (...) denotes an ensemble average, L is the total number of the random walk steps, xi is a single-step dis- 
placement and (s 2 ) is its mean-square value. In this derivation we neglected the off-diagonal term X)z^/' ( x l x l') or 
crossed correlations between successive single-step displacements in comparison with the diagonal term L (a; 2 ). This 
is because of the process definition, which invokes independent draw of steps. 

In Figure [31 the schematic illustration of the above considerations is shown by using a part of the trajectory or 
random walk realization consisting of hierarchically ordered steps (tq t 3 ,bo b 3 ), for j = 0, 1, 2. Herein, we neglected (i) 
fluctuation of the number of hierarchy levels js as well as (ii) their random succession. Thus we plotted the ordered 
trajectory within the time-space frame of coordinates. In fact, we made a transformation from the stochastic hierarchy 
to its deterministic counterpart. The explanation of the concept of black swans becomes now more convenient. 

We can easily derive the useful relation between the single-step mean-square displacement (x 2 ^ of the simulated 
trajectory and the maximal level, juAXi of the hierarchy contained in it. Indeed, the level juAX defines the extreme 
event or black swan by the pair of components (to t 3max , i>o b jMAX ). This quantity, for a large number of steps L>1 
or 2m AX > 1, is given by 

(x 2 ) « (b ) 2 (Mi^ (b 2 f + ( & 2)i + NiMAx-* (&2)2 + ... + J£ (pZ)mAx^ (39) 



where now (. . .) means an average over the random walk steps. That is, we used here a kind of an ergodic hypothesis. 
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FIG. 3: Schematic trajectory of hierarchically ordered steps (r 3 , b 3 ) presented, for simplicity, for N = 3, jmax = 2 and 
calibration parameters To = 1, &o = 1- Obviously, in our calculations we assumed jmax 2> 1. Herein, the extreme event or 
black swan is defined by the pair of components (t 3max=2 , yiMAx—^ directed by Jmax = 2. 



From Eqs. J3SJ and ([39]) we obtain 



, 01 „ (b 2 /N) jMAX+1 - 1 , ( i-4rnz(b 2 ) JMAX , for b 2 /N > 1 

(L) « (6 ) 2 iV— ( / J /N _ 1 « (M 2 | ^nLa*; for ^ < ! (40) 



(i-i/jv) 2//3 £2//} fo r < 2 



as well as 



(X 2 ) (L) « (fe„) 2 i-i/^ 1 ' ^ (41) 
I T^Tlv ^ for /3 > 2 



where the number of steps 



L w N JMAX + N^ MAX ~ l 4- i + /v° w atjmax C43\ 

1-1/N K J 

and the single-step displacement of black swan is defined by 

%M=&o6 3 ' M ". (44) 

Herein, the marginal case of /? = 2 is not considered. Apparently, for /3 < 2, i.e. if (^ 2 ) scales with L according to 
some power-law, (X 2 } is fully determined by x\ IAX . This is exactly what we need for illustration of our considerations. 
That is, the quantities (herein, the mean-square displacement of the process) which characterize the system are mainly 
expressed by the corresponding ones (herein, a single-step displacement) which define the black swans. 

Note that in the case of the WM-CTRW we have results analogous to those presented by Equations (|4"Tj) and (|4"2"j). 
However, their derivation is much more complicated in this case. The threshold property exhibited by these equations 
is typical of the behavior of other key quantities, like volatilities or correlation functions. 
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As expected, for Brownian motions (that is for the case of /3 > 2 in Equations (|40|) . (j4Tj) and (|42))). the factor 
preceding i in the second formula in (|41[) equals, in fact, to the corresponding factor present in Formula ([5]). Here, 
the absence of factor 2 is only caused by the absence of fluctuations of interevent times. 

Now, we are ready to answer the question concerning the distribution of the single-step displacements | x \ — &o W 
of the walker. This answer is based on the change of variables from j to | x |. Hence, the corresponding (normalized) 
distribution w(\ x |) takes the Pareto form 



where | x |> bo. In fact, Equation (|45]l holds only for | x |^> &o- It is well known [121 . 121l | that the power-law 
distribution of events leads to the Frechet distribution of extreme events. For asymptotic values of the argument, this 
distribution preserves the power-law of exponent f3 + 1 with the power-law correction to scaling of exponent f3. Note 
that Equation (|4"5)) is valid for all values of /3. That is, black swans are always present in the system. However, only 
for (3 < 2 their influence dominates. 

To conclude this discussion, we can say that the linear size of the walk (herein, the mean-square displacement) is 
controlled by extreme events or black swans if this size scales with the number of steps according to some power-law, 
i.e. if the walk is a kind of fractional random walk. Otherwise, the Brownian motions dominate and no influence of 
black swans is observed (although they are present in the system). That way the threshold functions as a discriminator 
of black swans. 

It is evident now, that we can identify an event as a sustained dragon-king if jd index of this event distinctly exceeds 
that of 3m ax ■ 



B. Comparison of theoretical predictions with simulations 



Herein, we restrict our simulations to an important example of a basic random walk given by a fractional Brownian 
motion. That is, we study a confined region of a superdiffusion phase defined by 1//3 only slightly smaller than 1/a. 
Note that inequality 1/(3 < 1/a is equivalent to v < 1, i.e. it corresponds to the case where the velocity of the walker 
for the higher hierarchy level is smaller. In this case, each moment of the arbitrary non-negative order of the (multi- 
step) displacement is finite for finite times (l7l . fisj ) . This choice of such moments arises from empirical evidences that 
these moments are always finite. Other choices, concerning other diffusion phases, would also be worth studying. The 
dragon-king is located "manually" inside the simulated time series in such a way that inequality At max *C tji is 
obeyed. 



1. Results for sustained dragon-king 

Now, we discuss the case where the multiplicative factor preceding in Equation (|26p is positive, which is easy to 

fulfil. That is, we consider inequality < 7^ ^1 — i-At/t t I ) ' This inequality is only slightly stronger than < 1, 

yet needed for the derivation of any of our expression for Cd(At) in the case of the appearance of the sustained 
dragon-king. 

In cumulative Table U we present data (in a form of four inverted pyramids of numbers), which define unnormalized 
statistics, S(j), of hierarchy levels js (cf. Section IY Aft . These data were obtained for fixed common parameters 
To = 1.0, t = 2.520, vq — 1.0, v — 0.992, and N = 4. They were used to prepare four trajectories within the WM- 
CTRW formalism in the presence of sustained dragon-kings. For instance, the row numbered by level j = 3 gives, at 
intersection with the second column, the number which says how many times (herein, it is 809231) this level appeared 
in the first trajectory. This trajectory contains the only sustained dragon-king defined by index jd = 13. This index 
is shown in Table U by bold number 1, at the intersection of row numbered by level j = 13 with the second column, 
again. The successive columns from three to five contain analogous unnormalized statistics but for an increasing jd 
values 6 , i.e. jd = 15, 17, and 19; that is, for increasing sustained dragon-kings. The common bottom of all hierarchies 
(represented in Table U by inverted pyramides of numbers) is placed at the level j ~ Jmax = 12, that is above any 
13 < j = j d < 19. 



Equations 1351 and 1361 precisely define the role of index used here. 
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TABLE I: Four unnormalized statistics S(j) of hierarchy levels js obtained for four sustained dragon-kings. 



Level j 


S(j) for j d = 13 


S(j) for j d = 15 


S(j) for h = 17 


S(j) for j d = 19 





51763445 


51439530 


49410801 


36182538 


1 


12948042 


12866869 


12360063 


9049915 


2 


3234819 


3214452 


3087567 


2260332 


3 


809231 


804047 


772246 


565401 


4 


202591 


201289 


193303 


141499 


5 


50583 


50521 


48326 


35374 


6 


12773 


12704 


12211 


8895 


7 


3162 


3141 


3012 


2212 


8 


811 


801 


765 


565 


9 


192 


191 


181 


128 


10 


48 


47 


43 


29 


11 


6 


6 


6 


4 


12 


5 


5 


4 


3 


13 


1 
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15 





1 








16 














17 








1 





18 














19 
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So, all sustained dragon-kings are marked in Table [J by the bold number 1. They are placed from the second to 
fifth column at intersections with the corresponding rows indexed by levels from j = jd = 13 to j = jd = 19. That is, 
these levels systematically move away from the common bottom of hierarchies. 

Importantly, the left part of trajectory (preceding the dragon-king appearance) is common for all dragon-kings. 
The right border of this part is fixed defining the beginning of a dragon-king. Because the total duration time, ttot, 
of all trajectories is the same, the duration time, tn, of the part of trajectory placed on the right-hand side of the 
dragon-king decreases as the duration time td increases (i.e. when the index jd increases). Therefore, the statistics 
S{j) of hierarchy levels js, shown in Table HI decreases as the level defining the sustained dragon-king, jd, is rised 
(because then the corresponding trajectory or time series, placed on the right-hand side of the dragon- king, is shorter). 
Therefore, for example, the number 128 placed at intersection of the fifth column and the row denoted by the index 
level j = 9 is distinctly smaller than the number 192 placed at the same row but at the intersection with the second 
column. We hope that Table U well illustrates the hierarchical structure of any (long) trajectory simulated within 
the WM-CTRW formalism. Moreover, the corresponding localisations of the sustained dragon-kings in the space of 
hierarchy levels relative to the inverted pyramids are also well delineated. 

In Figure[H we compare the prediction of Formula ([2U]) (thin solid curves) with results of simulation (dispersed thick 
solid curves) for four different values of td, namely, td/ AtMAX = 1-653, 10.496, 66.651, and 423.263, which correspond 
to jd = 13, 15, 17, and 19, respectively. Note that the current width of the time-window, At, (called also time lag) 
ranges from At — dt up to At = AtMAX = 10 5 dt with the time step, dt — 1, while t to t = 1400 At max is the same 
for all statistics S{j) shown in Table [I] The sustained dragon-kings' time lags were calculated for the same value of 
t = 2.520 and the single-step displacements were calculated for the common b = vt = 2.50 (where bo = vqTq = 1.0). 

The upward convexity of the curves in Figure 0] is due to the presence of the corresponding sustained dragon-king. 
As we expected, the agreement shown in Figure 0] between the prediction of Formula ([2"?J]) and data obtained from 
simulations becomes better the further the dragon-king is located from the top of the hierarchy (see the location of 
bold number 1 in Table []}. The best agreement is obtained for the largest jd = 19. In other words, the dragon-kings 
defined by jd smaller than 19 slightly positively deviate their corresponding velocity autocorrelation functions, Cd{At), 
from their simulational counterparts. 

The quantity C(At), controlled by black swans and given by Expression ([TT]) . is also plotted in Figure H] (the dashed 
curve) as a reference VAF, i.e. this VAF was calculated for the absence of a dragon-king. 
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FIG. 4: Comparison of the prediction of Equation (|26|l (solid thin curves) with results of simulations (dispersed solid thick 
curves) for four different values of td/ AtMAX = 1.653, 10.496, 66.651, and 423.263 which correspond to jd = 13, 15, 17, and 19, 
respectively. The dashed curve represents the prediction of Equation (|11|) . i.e. the prediction for the time series in the absence 
of a dragon-king. All curves were calculated for the same values of r = 2.52, v — 0.992, and TV = 4. Notably, theoretical 
predictions for jd — 17 and 19 are almost indistinguishable (within the resolution of the plot) from results of the corresponding 
simulations for the whole range of time lag At. The presence of the sustained dragon-king twists upward the curve deviating 
it from the straight line (in the log — log plot). 



2. Results for the shock dragon-king 

In Figure^ we present results of simulations of the VAF for three different values of the shock size Xd = 0.41, 2.44 
and 5.26 [xlO 6 ] which corresponds to X,i/xmax = 1-90, 11.3, and 24.4, respectively, where Xmax is the maximal 
spatial value of the random walk's single step belonging to the simulated hierarchical WM-CTRW trajectory in the 
absence of a dragon-king. Remarkably, all these shocks are fully exogenous as they were taken from outside of 
the spatio-temporal structure of time series or random walks. The simulated trajectories have also the total time 
t to t = 1400 AtMAX and, except for the presence of the dragon-king, all trajectories are identical. 

A striking property of simulated VAFs is their dispersion behaving like a certain instability. Therefore, it is 
more convenient to use a formula that only describes the dispersion of the data. In principle, such a formula could 
be obtained by replacing the sum of velocities vt, + vr in Formula (|32p by its dispersion a = yJ (Jvl + v r) 2 ) = 
V / 2a/°'v + C(2At), where a% = (i>f,) = (vr)- This is allowed because we assumed, for simplicity, that a shock does 
not change the type of the random walk. 

Moreover, our approach allows us to study a more realistic case, e.g. in which dispersion is smaller than a defined 
above. Hence, we propose a more flexible stationary formula 

Ci (A i )= C (A t)± ^,^-(^) 2 , (46) 
where 07 d =' \J fa% + C(2At) and the phenomenological factor or weight, < / < 1, is the same for all trajectories. 
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FIG. 5: Comparison of the prediction of Formula (|46|) (dotted and dashed-dotted curves as well as grey solid one) with results 
of simulations (corresponding regions having different greyness, additionally marked by a, b and c) for three different values of 
Xd = 5.26, 2.44, and 0.41 [xlO 6 ]. The black solid curve shows the prediction of Formula (|lip . i.e. the prediction for the time 
series simulated in the absence of a shock dragon-king; the corresponding simulated VAF is given for this time series by the 
innermost (dark) region. All curves were obtained for values of r = 2.52, v = 0.992, and N = 4, the same as for the case of 
the sustained dragon-king. 



As it is seen, also in the case of the shock dragon-king we transformed the non-stationary Expression (J35]) to more 
useful stationary Expression (|4"?J]) . 

Comparison of the prediction of Equation (|46p with the data obtained by simulation is shown in Figure [5] In this 
figure, only small deviations are seen for the choice of the factor / = 0.30. The data scatter is reasonably small but 
it increases with the increase of the ratio Xd/ttot- The origin of this scatter comes from fluctuations of the simulated 
trajectory, unfortunately resulting also in a spontaneous artificial trend (e.g. as a deviation from the power-law in the 
absence of the dragon-king). Additionally, this trend can be supported by the finite size of the simulated time series. 

Note that further simplification of Equation ([4"fJ)) is also possible 

C d (M) = C(At) ±V2a f ^~ (^) , (47) 

Hot \Hot J 

if the strong but reasonable inequality At max /ttot <C 1 is obeyed. 



VI. SUMMARY AND CONCLUDING REMARKS 



In the present work we discussed the following issues. 
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(i) We considered the influence of two distinctive types of dragon-kings on the velocity autocorrelation function, 
namely the sustained dragon-king and shock dragon-king. 

(ii) By simulations and by theoretical analysis, we found that the dragon-king influence decisively changes the 
original VAFs calculated, for instance, within the hierachical Weierstrass-Mandelbrot Continuous-Time Random 
Walk formalism for a wide range of time intervals. The influence of both types of dragon-kings is well pronounced 
but quite different (cf. the corresponding plots shown in Figures @] and [5]) . Remarkably, the results obtained by 
simulations agree well with the corresponding predictions of our simple theoretical Formulas and P5|) (see 
again Figures 2] and [5]) . 

(hi) Furthermore, several intermediate formulas, e.g. (|20|) . (|22)) . or (|30|) . derived in this work for Cd(At), can be 
applied to more complex cases where, for instance, (a) after the dragon-king appearance the random walk is 
changed, (b) the random walk with drift is considered and (c) dragon-kings cluster in the system. 

As it is apparent from Figure 0J the presence of the sustained dragon- king throws the system away from the state 
controlled by black swans. That is, the 'finger-print' of this dragon-king is significant because it convexly twists 
upward the autocorrelation curve dependence on the time lag and distinctly deviates it positively from the power-law 
generated by black swans (cf. Section IV A 21 for details). This deviation is one of the main features which distinguishes 
Cd(At) from usual C(At). This difference allows one to distinguish a power-law relaxation, controlled by black swans, 
from the decay controlled by the sustained dragon-king. Moreover, this difference can be so distinct that it can 
effectively be applied as a tool to detect the sustained dragon-king from an empirical time series. 

The plot in Figure [3] shows that the shock dragon-king also significantly changes C(At). However, this change is 
different than the change caused by the sustained dragon-king. This change is also well seen by direct comparison of 
Formulas (|2"6"|) and (|4"u]) . Noticeably, the scatter of the data shown in Fig. [S] indicates some instability of the system 
after the appearance of the shock dragon-king. Indeed, if the empirical VAF reveals such an instability then we can 
anticipate that the corresponding empirical time series contains the shock dragon-king. That is, such an anomalous 
VAF indicates shocks. 
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